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The theoretical status of the 'proton spin' effect is reviewed. The conventional QCD parton model analysis 
of polarised DIS is compared with a complementary approach, the composite operator propagator-vertex (CPV) 
method, each of which provides its own insight into the origin of the observed suppression in the first moment of 
g\ . The current status of both experiment and non-perturbative calculations is summarised. The future role of 
semi-inclusive DIS experiments, in both the current and target fragmentation regions, is described. 



1. Introduction 

It is now nearly a decade since the EMC collab- 
oration announced measurements H of the first 
moment of the polarised proton structure func- 
tion g\ which in a simple parton picture could 
be interpreted as showing that the quarks carry 
only a small fraction of the spin of the proton. 
This 'proton spin' problem has been the focus 
ever since of an extraordinary experimental and 
theoretical effort to extend and improve the data 
and provide a thorough understanding of the phe- 
nomenon in terms of QCD. 

Perhaps the 1998 Montpellier conference will 
mark the tenth anniversary of the EMC paper 
with a comprehensive review of all this work. In 
this ninth anniversary review, I will instead take 
the opportunity to present a more individual look 
at what the 'proton spin' effect means, how it 
can be understood in non-perturbative QCD, and 
what future experiments it suggests. 

The central theme of this review is the com- 
parision of two complementary approaches to the 
description of deep inelastic scattering (DIS)- 
the parton model and the composite operator 
propagator- vertex (CPV) method J2|-^|]. We will 
show how each provides its own insight into the 
'proton spin' effect - the first from a quark-gluon 
constituent point of view and the second combin- 
ing a holistic description of the proton with non- 
perturbative, target-independent, QCD physics. 



The status of non-perturbative calculations us- 
ing QCD spectral sum rules (QSSR) and lattice 
methods is then briefly reviewed, followed by a 
discussion of the current experimental situation 
for gf(x, Q 2 ) and the important unresolved ques- 
tions concerning the small x region. Finally, the 
present and future role of semi-inclusive DIS ex- 
periments, both in the current and target frag- 
mentation regions, is considered. A recent pro- 
posal^] for testing the target-independent sup- 
pression mechanism suggested by the CPV ap- 
proach is described and predictions for cross sec- 
tion moment ratios for eN — > ehX are given. 

2. The sum rule for 1^ and the OZI rule 

The polarised structure functions are mea- 
sured in polarised DIS (see Fig. 1), either /ip — > 
LiX (EMC,SMC) or ep -> eX (SLAC, HERMES). 
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Fig.l Inclusive polarised DIS scattering. 

g\ is extracted from the polarisation asymmetry 
of the cross section according to: 
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where Q 2 



-q 2 and x = Q 2 /2p2-q are the 



Bjorken variables, y = Q 2 jxs and Yp = (2 — y)/y. 
(The dots denote terms of 0{M 2 x 2 /Q 2 ), includ- 
ing the second polarised structure function g\ for 
which interesting new measurements are becom- 
ing available ||.) 

On the theoretical side, g\ is determined by 
the proton matrix element of two electromagnetic 
currents carrying large spacelike momentum. The 
sum rule for the first moment Y\ is derived using 
the twist 2, spin 1 terms in the operator product 
expansion (OPE) for the currents: 



V3 " 



;C!(a s )Al (2) 



where the Wilson coefficients C^ s (a s ) and 
Cf (a s ) are now both known to 0(a 3 ) M. It reads: 
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C?a°(Q 2 ) (3) 



Here, a 3 , a 8 and a°(Q 2 ) are the form factors in 
the forward proton matrix elements of the renor- 
malised axial current, i.e. 



(p,«|A£|p,«) =8, 
(p,s\Al\p,s) =8, 
( P ,s\Al\p 7S )= s y(Q 2 ) 
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where p^ and s M are the momentum and polari- 
sation vector of the proton. 

Because of the chiral Ua(1) anomaly, d^A^ — 
2nfQ ~ 0, the flavour singlet current A^ is not 
conserved. It is renormalised and mixes with the 
topological density. Defining the bare operators 
Kb = E?7m759 and Q B = ^e^trG^G^, 
we have (for nt flavours) 



Al = ZA% 



Q = Qb 



2n, 



(1 - Z)d»Al 



(5) 



where Z is a divergent renormalisation constant. 
The associated anomalous dimension 7 is known 
to 3 loops Q. Matrix elements of A ^ therefore 



have a non-trivial renormalisation group (RG) 
scale dependence governed by 7. In particular, 



(6) 



where t = lnQ 2 /A 2 , so that the singlet axial 
charge a°(Q 2 ) is scale dependent. This is cru- 
cial to understanding the 'proton spin' effect in 
QCD. 

The axial charges a 3 and a 8 are known in terms 
of the F and D constants (a 3 = F + D, a 8 = 
3F — D) found from neutron and hyperon beta 
decays. The interest of the sum rule centres on 
the flavour singlet axial charge a°(Q 2 ). In the 
absence of an alternative experimental derivation 
of a°(Q 2 ), the simplest ansatz is to assume that 
it obeys the OZI (Zweig) rule^, i.e. a°(Q 2 ) = a 8 . 
This gives a theoretical prediction for which is 
known as the Ellis- Jaffe sum rulc(l^]. This is now 
known to be violated (see sect. 6), with a°(Q 2 ) 
strongly suppressed relative to a 8 . This is what 
is known as the 'proton spin' problem. 

In fact, it is not at all surprising that the OZI 
rule should fail in this casejll]]. The first clue is 
the anomaly-induced scale dependence of a°(Q 2 ). 
If the OZI rule were to hold, at what scale should 
it be applied? For the same reason, it is imme- 
diately clear that a°(Q 2 ) cannot really measure 
spin. Moreover, it is known that the pseudovec- 
tor and pseudoscalar channels are linked through 
the Goldberger-Treiman relations^). Since large 
anomaly-induced OZI violations are known to be 
present in the pseudoscalar channel (L0t(l) prob- 
lem, 7/ mass, etc.) it is natural to find them also 
for a°(Q 2 ) in the pseudovector channel. 

While this immediately resolves the 'proton 
spin' problem, clearly we want to understand the 
origin of the suppression in a°(Q 2 ) much more 
deeply. The following two sections describe two 
complementary approaches to this question - the 



3 The OZI limit of QCD is defined (see ref . ) as the trun- 
cation of full QCD in which non-planar and quark-loop 
diagrams are retained, but diagrams in which the exter- 
nal currents are attached to distinct quark loops (so that 
there are purely gluonic intermediate states) are omitted. 
This is a more accurate approximation to full QCD than 
either the leading large l/N c limit, the quenched approx- 
imation (small nj at fixed JV C ) or the leading topological 
expansion (N c — > 00 at fixed rif/N c ). 
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conventional QCD parton_model and the CPV 
method developed in refs. 

3. The QCD parton model 

In the most simple parton model, where the 
proton structure for large Q 2 DIS is described by 
parton distributions corresponding to free quarks 
only, the polarised structure function is given by 



9 P M = 



1 nf 



Aqi(x) 



(7) 



where Aq t (x) = qf (x) + <j+ (x) - qT {x)-q~ (x) is 
the difference of the distributions of quarks with 
helicities parallel and antiparallel to the nucleon 
spin. It is convenient to work with the flavour 
non-singlet and singlet combinations: 



Aq NS {x) 
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Aq S (x) =J2 A <H(x) 



(8) 



In this model, the first moment of the flavour sin- 
glet quark distribution Aq s = L dx Aq s (x) can 
indeed be identified as the sum of the helicities 
of the quarks. Interpreting the structure function 
data in this model then leads to the conclusion 
that the quarks carry only a small fraction of the 
spin of the proton - the 'proton spin' problem. 

However, this simple model leaves out many 
important features of QCD - gluons, RG scale de- 
pendence, the chiral Ua(X) anomaly, etc. When 
these effects are included, in the QCD parton 
model, the naive identification of Aq s with spin 
no longer holds and the experimental results for 
g\ are readily accommodated. 




Fig. 2 QCD parton model interpretation of DIS. The single 
lines are partons, which may be quarks or gluons. 



The QCD parton model picture of DIS is shown 
in Fig. 2. The polarised structure function is writ- 
ten in terms of both quark and gluon distributions 
as follows: 



+ C S QAq S (y,t) 



NS 



(*)Aq NS (y,t) 



C[-)Ag(y,t) (9) 



where C s , C 9 and C NS are perturbatively calcu- 
lable functions related to the Wilson coefficients 
in sect. 2 and the quark and gluon distributions 
have a priori a t = lnQ 2 /A 2 dependence. 

The RG evolution (DGLAP) equations for 
these polarised distributions are: 



Aq NS (x,t) 



dt 



2tt 



dy 

y 
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Aq NS (y,t)(10) 



and, abbreviating the notation of (10), 
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showing the mixing between the singlet quark and 
the gluon distributions. The splitting functions 
P are also calculable in perturbative QCD, their 
moments being related to the anomalous dimen- 
sions of the series of increasing spin operators ap- 
pearing in the OPE (2). 

In this language, the first moment sum rule for 
gi reads: 



r?(Q 2 



C? S Aq NS 



CfAg 5 + Cf Ag (12) 



where Aq NS , Aq s and Ag are the first moments 
of the above distributions. Comparing with (3), 
we see that the axial charge a°(Q 2 ) is identified 
with a linear combination of the first moments of 
the singlet quark and gluon distributions. It is 
often, though not always, the case that the mo- 
ments of parton distributions can be identified in 
one-to-one correspondence with the matrix ele- 
ments of local operators. The polarised first mo- 
ments are special in that two parton distributions 
correspond to the same local operator. This adds 
an extra subtlety to the identification. 

The RG equations for the first moments of 
the parton distributions follow immediately from 
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(10,11) and depend on the matrix of anomalous 
dimensions for the lowest spin, twist 2 opera- 
tors. This introduces a renormalisation scheme 
ambiguity. The issue of scheme dependence has 
been studied thoroughly by Ball, Forte and Ri- 
dolfi |]12| a nd an excellent summary can be found 
in ref.[|l3|. It is shown there that it is possible to 
choose a scheme known as the Adlcr-Bardeen or 
AB scheme (strictly, a class of schemes |l2|,[l3| ) for 
which the parton distributions satisfy the follow- 
ing RG equations: 



d a s 
~dt 2?r 



(13) 



with the implication Cf = —Uf^-Cf. It is then 
possible to make the following identifications with 
the axial charges: 

a 3 = Au - Ad 



= Au + Ad - 2As 



a°(Q 2 ) = Aq s -n f ^-Ag(Q 2 



(14) 



where Au = <ix(Au(x, t) + Au(x, i)) etc. No- 
tice that in the AB scheme, the singlet quark dis- 
tribution Aq s (which is often written as AE) is 
scale independent. All the scale dependence of 
the axial charge a°(<5 2 ) is assigned to the gluon 
distribution Ag(Q 2 ). 

This was the identification originally intro- 
duced for the first moments by Altarelli and 
Ross|l4|]. We emphasise that (13) is true only 
in a particular renormalisation scheme (the AB 
scheme) and that it is only in this scheme that 
the identifications (14) hold. 

In this picture, the Ellis- Jaffe sum rule follows 
from the assumption that in the proton both As 
and Ag(Q 2 ) are zero. This is equivalent to the 
naive OZI approximation a°(Q 2 ) = a 8 described 
above. Clearly, given the RG scale dependence of 
a°(Q 2 ), this assumption is in contradiction with 
QCD where the anomaly requires a°(Q 2 ) to scale 
with the anomalous dimension 7. 

Since neither AS nor Ag(Q 2 ) are currently 
measurable in other processes, the parton model 



is unable to make a quantitative prediction for the 
first moment T^. While the model can accomo- 
date the observed suppression, it cannot predict 
it. 

An interesting conjecture, proposed in the orig- 
inal paper of Altarelli and Ross 14 . is that the 
observed suppression in a°(Q 2 ) is due overwhelm- 
ingly to the gluon distribution Ag(Q 2 ). If so, the 
strange quark distribution As ~ in the proton 
and so AS ~ a 8 . This is entirely plausible be- 
cause it is the anomaly (which is due to the glu- 
ons and is responsible for OZI violations) which 
is responsible for the scale dependence in a°(Q 2 ) 
and Ag(Q 2 ) whereas (in the AB scheme) AS is 
scale invariant. The essence of this conjecture will 
reappear in the next section where we describe 
the CPV method. 

To test this conjecture, we need to find a way to 
measure Ag(Q 2 ) itself, rather than the combina- 
tion a°(Q 2 ). One possibility (see also sect. 6) is to 
exploit the different scaling behaviours of Aq s (x) 
and Ag(x 1 Q 2 ) to distinguish their contributions 
in measurements of g\(x, Q 2 ) at different values 
of Q 2 . A second is to extract Ag(x, Q 2 ) from pro- 
cesses such as open charm production, 7*3 — > cc, 
which will be studied in various forthcoming ex- 
periments at COMPASS, RHIC, etc. 



4. The CPV method for DIS 

This approach to sum rules in DIS was devel- 
oped in collaboration with S. Narison and 
G. Veneziano in a series of papers on the 'pro- 
ton spin' effect. The starting point, as described 
above, is the use of the OPE in the proton matrix 
element of two currents. This gives the standard 
form for a generic structure function moment: 

[ dx x"-^(x; Q 2 ) = V C"(Q 2 )(p|Or(0)b)(15) 
Jo 



where Of are the set of lowest twist, spin n oper- 
ators in the OPE and C™(Q 2 ) the corresponding 
Wilson coefficients. In the CPV approach, we 
now factorise the matrix element into the prod- 
uct of composite operator propagators and vertex 
functions, as illustrated in Fig. 3. 
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Fig. 3 CPV description of DIS. The double line denotes 
the composite operator propogator and the lower blob the 
'1PP vertex. 

To do this, we first select a set of composite op- 
erators Oi appropriate to the physical situation 
and define vertices 1^ as '1PF with respect to 
this set. Technically, this is achieved by intro- 
ducing sources for these operators in the QCD 
generating functional, then performing a Legen- 
dre transform to obtain an effective action r[Oj]. 
The 1PI vertices are the functional derivatives of 
r[Cy. The generic structure function sum rule 
(15) then takes the form 



dx x n ~ x F{x,Q 2 ) = 



This decomposition splits the structure func- 
tion into three pieces - first, the Wilson coef- 
ficients Cj n \Q 2 ) which control the Q 2 depen- 
dence and can be calculated in perturbative QCD; 
second, non-perturbative but target- independent 
QCD correlation functions (composite operator 

propagators) (0\T O^ 6i\Q); and third, non- 
perturbative, target-dependent vertex functions 
^Oipp describing the coupling of the target proton 
to the composite operators of interest. The ver- 
tex functions cannot be calculated directly from 
first principles. They encode the information on 
the nature of the proton state and play an analo- 
gous role to the parton distributions in the more 
conventional parton picture. 

As emphasised in refs.|§-Q|, it is important to 
recognise that this decomposition of the matrix 
elements into products of propagators and proper 
vertices is exact, independent of the choice of the 
set of operators Oi. In particular, it is not neces- 
sary for Oi to be in any sense a complete set. All 



that happens if a different choice is made is that 
the vertices ^o iPP themselves change, becoming 
'1PI' with respect to a different set of composite 
fields. Of course, while any set of Oi may be cho- 
sen, some will be more convenient than others. 
Clearly, the set of operators should be as small as 
possible while still capturing the essential physics 
(i.e. they should encompass the relevant degrees 
of freedom) and indeed a good choice can result in 
vertices T^, pp which are both RG invariant and 
closely related to low energy physical couplings, 
such as g-xNN or (7 7r77 [|). In this case, (16) pro- 
vides a rigorous relation between high Q 2 DIS and 
low-energy meson-nucleon scattering. 

For the first moment sum rule for g\, it is most 
convenient to use the chiral anomaly immediately 
to re-express a°(Q 2 ) in terms of the forward ma- 
trix element of the topological density Q, i.e. 



a°(Q 2 ) 



1 

2M 



2n s {p\Q\ P ) 



(17) 



Our set of operators Oi is then chosen to be 
the renormalised flavour singlet pseudoscalars Q 
and $5 where, up to a crucial normalisation fac- 
tor, the corresponding bare operator is $55 = 

9759- This normalisation factor is chosen such 
that in the absence of the anomaly, or more pre- 
cisely in the OZI limit of QCD (see footnote 3), 
$5 would have the correct normalisation to cou- 
ple with unit decay constant to the U(l) Gold- 
stone boson which would exist in this limit. This 
also ensures that the vertex is RG scale indepen- 
dent 01. We then have 



1 1 

1 " l " 3 ~ 9 2M" 



r? „,„„ = l^n/Cf (a.) x 



\TQ Q\0)r Qpp + (0\TQ $ 5 |o)r$ 5j 



(18) 



where the propagators are at zero momentum and 
the vertices are 1PI wrt Q and $5 only. This is 
illustrated in Fig. 4. 



XQ 



XQ 



Fig. 4 CPV decomposition of the matrix element (p|Q|p). 
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The composite operator propagator in the first 
term is the zero-momentum limit of the QCD 
topological susceptibility x(k 2 ), viz. 



X (k 2 ) = / dxe lk - x i(0\T Q(x) Q(0)|0) 



(19) 



The anomalous chiral Ward identities show that 
X(0) vanishes for QCD with massless quarks, in 
contrast to pure Yang-Mills theory where x(0) is 
non-zero. Furthermore, it can be shownQ that 
the propagator (0\T Q $5|0) at zero momentum 
is simply the square root of the first moment of 
the topological susceptibility. We therefore find: 

rf sing = l^ 2n f c f( a s) VS r* 5PP (20) 



The quantity y'x'(O) is not RG invariant and 
scales with the anomalous dimension 7, i.e. 



dt 

On the other hand, the proper vertex has been 
chosen specifically so as to be RG invariant. The 
renormalisation group properties of this decom- 
position are crucial to the resolution proposed in 
ref.[g|| of the 'proton spin' problem. 

The proposal of ref. ^,|) is that we should ex- 
pect the source of OZI violations to lie in the RG 
non-invariant, and therefore anomaly-sensitive, 
terms, i.e. in x'(0)- Notice that we are us- 
ing RG non-invariance, i.e. dependence on the 
anomalous dimension 7, merely as an indicator 
of which quantities arc sensitive to the anomaly 
and therefore likely to show OZI violations. Since 
the anomalous suppression in is assigned to 
the composite operator propagator rather than 
the proper vertex, the suppression is a target in- 
dependent property of QCD related to the chi- 
ral anomaly, not a special property of the proton 
structure^. 

To convert this into a quantitative prediction 
we use the OZI approximation for the vertex 

4 Other models whichjxeedict target independence exist in 
the literature. In ref. [|15| it is suggested that the suppres- 
sion in a°(Q 2 ) may be due directly to non-pertuxbative 
effects in 7 at low scales. In the model of ref. pq], the 
axial charge is related to the nucleon couplings of the 
pseudovector mesons, which are relatively uncertain ex- 
perimentally, rather than using the Goldberger-Treiman 
relations to compare with the pseudoscalar sector. 



r$ 5 pp and a QCD spectral sum rule estimate of 
the first moment of the topological susceptibility. 
We find||, for n f = 3, 



VxW) 



Q 2 = 10GeV 2 



= 23.2 ± 2.4 MeV 



(22) 



This is a suppression of approximately a factor 
0.6 relative to the OZI value f^/y/6. 

Our final result, in the chiral limit, is then 



a°(Q 2 = lOGeV 2 ) = 0.35 ±0.05 
from which we deduce 



Q2 = ioGcV 2 



= 0.143 ±0.005 



(23) 



(24) 



This is to be compared with the Ellis-Jaffe (OZI) 
prediction of a = 0.58. 

The complementary nature of the QCD par- 
ton model and CPV methods is now clear. Both 
involve at present incalculable non-perturbative 
functions describing the proton state - the quark 
and gluon distributions in the parton picture and 
the 1PI vertices in the CPV method. Both ex- 
hibit a degree of universality - the same parton 
distributions may be used in different QCD pro- 
cesses such as DIS or hadron-hadron collisions, 
while the vertices (when they can be identified 
with low-energy couplings) also provide a link be- 
tween high Q 2 DIS and soft meson-nucleon inter- 
actions. 

One of the main advantages of the CPV 
method is that some non-perturbative informa- 
tion which is generic to QCD, i.e. independent 
of the target, is factored off into the compos- 
ite operator propagator. This allows us to dis- 
tinguish between non-perturbative mechanisms 
which are generic to all QCD processes and 
those which are specific to a particular target. 
As explained above, our contention is that the 
anomalous suppression in the first moment of g\ 
is of the first, target-independent, type. This 
target-independence conjecture can in principle 
be tested by DIS with non-nucleon targets. This 
may effectively be realised in semi-inclusive DIS 
(see sect. 7 and ref.[[|). 

Both the parton and CPV methods allow a nat- 
ural conjecture in which the origin of this sup- 
pression is attributed to 'glue' - either through 
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a large polarised gluon distribution Ag(Q 2 ) in 
the parton description or due to an anomalous 
suppression of the first moment of the topologi- 
cal susceptibility yY(0) in the CPV description. 
These conjectures are based on assumptions that 
the appropriate RG invariant quantities, AS or 
r$ 5j , p , obey the OZI rule. The motivation for this 
is particularly strong in the CPV case, where it is 
supported by a range of evidence from low-energy 
meson phenomenology in the Ua(1) channel. 

The two approaches therefore provide related, 
but complementary, insights into the nature of 
the 'proton spin' effect. Perhaps the most impor- 
tant message of this review is that both insights 
are needed and both methods have a full part to 
play in understanding this intriguing and subtle 
phenomenon. 

5. Non-perturbative results QSSR and 
lattice 

The challenge to non-perturbative calcula- 
tional methods in QCD is therefore to evalu- 
ate a°(Q 2 ), either directly through the matrix 
element or, using the anomaly, from 

(p\Q\p)- Alternatively, if we accept the conjec- 
ture proposed above, we can deduce a°(Q 2 ) from 
a calculation of x'(0). In any case, the topological 
susceptibility is a fundamental correlation func- 
tion in QCD of great importance in a variety of 
contexts and deserves to be studied in its own 
right. 

These non-perturbative calculations can be 
performed either using lattice gauge theory or 
the method of QCD spectral sum rules (QSSR). 
It must be stressed that to be meaningful, these 
calculations must be performed in full QCD, or 
at the very least in an approximation that in- 
corporates chiral Ua(1) breaking by the anomaly. 
The occurrence of massless singularities is quite 
different in truncations such as large N c or the 
quenched approximation (see footnote 3) where 
the anomaly is not properly included, due to the 
erroneous appearance of a £/a(1) Goldstone bo- 
son. 

Lattice calculations of a°(<5 2 ) from 
have been performed by several authors and are 
reviewed in ref.jnj. Since only quenched config- 



urations have been used, the results cannot give 
the true a°(Q 2 ). However, a major step towards 
the true answer can be made by including explic- 
itly OZI-violating 'disconnected' diagrams with 
purely gluonic intermediate states as well as the 
OZI respecting 'connected' diagrams (see Fig. 5). 




P P P P 

(a) (b) 

Fig. 5 'Connected' (a) and 'disconnected' (b) contributions 
to the matrix element 

Using this method, Dong, Lagae and Liu llSl 
[19) find a = 0.25 ± 0.12, and confirm that 
the contribution from the connected diagrams 
(a° conn = 0.62 ± 0.09) reproduces the OZI expec- 
tation. It is not clear, however, how much this re- 
sult will change if non-quenched gauge configura- 
tions including the effects of dynamical fermions 
are used. 

For (p\Q\p), so far only quenched calcula- 
tions have been performed |p7|, although the Pisa 
group (2^] are currently running simulations in full 
QCD. Existing calculations of proton matrix el- 
ements using QSSR have also not fully incorpo- 
rated OZI breaking. 

The first moment of the topological suscepti- 
bility, x'(0), has been evaluated using QSSR in 
ref.(|) and the result is quoted above in (22). This 
shows Va(0) to b e suppressed to around 0.6 of 
the OZI value. The method and calculational de- 
tails are described in full in ||. 

The application of QSSR to the Ua(1) sector 
of QCD has been criticised repeatedly by Ioffc 
(see e.g. (H)) on two grounds: (i) that there are 
important neglected contributions from 'instan- 
tons', i.e. higher dimensional condensates, and 
(ii) that when the strange quark mass is in- 
cluded, QSSR results for current-current correla- 
tors show quite unrealistic SU (3) breaking. Nei- 
ther criticism is valid, and both will be refuted 
in detail in a forthcoming paper(22). In fact, 
the stabilisation scale in the calculation of ref. || , 
viz. t _1 ~ 2GeV 2 , is sufficiently big for higher 
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dimensional condensates to be suppressed, and 
indeed the calculation does display the hierarchy 
of gluonium to light meson hadronic scales an- 
ticipated by ref.^3|. Ioffe's second criticism is 
based on a calculation where radiative correc- 
tions are not properly implemented, giving a false 
result. The correct results are given in ref. p2| , 
where we extend our previous analysis systemat- 
ically beyond the chiral limit using a new set of 
generalised Goldberger-Treiman relations. 

x'(0) is a particularly difficult correlation func- 
tion to calculate on the lattice, requiring algo- 
rithms that implement topologically non-trivial 
configurations in a sufficiently fast and efficient 
way. However, very preliminary results from the 
Pisa group |^] of calculations in full QCD with 
dynamical quarks indicate a value of y/ x'(0) ~ 
f 6 ± 3MeV. Given the preliminary nature of the 
lattice simulations, the agreement with the QSSR 
result (22) is encouraging and further results from 
the lattice are eagerly awaited. 

6. Experiment and the small x region 

The most recent published results from the 
SMC collaboration on g\ are given in ref. |p3 . For 
the first moment, SMC quote: 

ft = 0.136 ± 0.013 ± 0.009 ± 0.005(25) 

Q 2 = 10GcV 2 

where the last error is a theoretical one, related to 
the evolution of the measured data to a standard 
Q 2 . This implies 



a°(Q 2 = lOGeV 2 



0.28 ±0.16 



(26) 



This is to be compared with the prediction (23,24) 
obtained in sect. 4. 

Although this agreement is very promising and 
suggests strongly (see also below) that the ex- 
planation of the 'proton spin' effect in terms of 
a suppression in the topological susceptibility is 
correct, there is an important uncertainty in the 
presentation of this data. In fact, SMC only 
take measurements in the region x > 0.003. The 
contribution of the small x region to T\ is esti- 
mated in ref. [^i) by a dubious Regge extrapola- 
tion, which gives only a very small addition of 
0.0042 ± 0.0016. Ball, Forte and Ridolfipl (see 



also J25|]) have proposed an alternative small x 
extrapolation, using the data to fit the parton 
densities at small a; at a low Q 2 scale, then using 
the perturbative QCD evolution equations to de- 
duce the form of g\{x, Q 2 — lOGeV 2 ). Examples 
of their fits are shown in Fig. 6. 




solid: fit A 
dashed: fit B 
dotted: fit C 
dotdashed: fit D 



Fig.6 QCD evolved g\ (x, Q 2 = lOGeV 2 ) at small X. 

The important feature of the evolution equa- 
tions is that while AS(i, Q 2 ) falls with increasing 
Q 2 , the polarised gluon distribution Ag(x,Q 2 ) 
rises. The net effect is that gP(x,Q 2 ) is driven 
strongly negative at Q 2 = lOGeV 2 for sufficiently 
small x. This gives a potentially large negative 
contribution to the first moment T^, but with rel- 
atively large errors. 



This summer, SMC [26 271] have released new, 
still preliminary, data on g\ based on the 1996 
run. See Fig. 7. 



gj at measured Q 



preliminary 



A 93+96 combined 
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10 



Fig.7 Preliminary SMC data for g x 



9 



It shows two significant features compared with 
the older data - the smallest x point has dropped 
appreciably, while there is a small rise in the 
medium x data points. 

xgf 
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Fig. 8 Preliminary SMC data for xg\ including 'Regge' 
and 'QCD' small X extrapolations. 

In presenting these results (see Fig. 8), SMC 
have given fits to the small x region using both the 
'Regge' and 'QCD' extrapolations. They quote 
the following results. For the measured range 
(including an uncontroversial x > 0.7 extrapo- 
lation) : 



f dx gft(x;Q 2 = lOGeV 2 ) 

^0.003 



0.146 ±0.006 ±0.009 ± 0.005 (27) 
while with small x extrapolations: 



r? 



Regge 



= 0.149 ±0.006 ±0.009 ±0.005(28) 

(29) 



Q 2 = 10GcV 

a°(Q 2 = 10GeV 2 ) i?e " e = 0.41 ± 0.11 



and 



r? 



QCD 



0.135 ± 0.006 ± 0.009 ± 0.011(30) 



Q 2 = 10GeV 2 

°(Q 2 = 10GeV 2 ) QCD = 0.27 ± 0.15 



(31) 



Notice that these results are for g\ taken from 
proton data alone. The values of T\ and a for gf 
taken from combined proton and deuteron data 
are systematically lower p6| , p7[ . 



The small x region is therefore an impor- 
tant challenge to experimentalists. Resolving the 
above uncertainty in and a will enable a rigor- 
ous test of the CPV conjecture and the link with 
the topological susceptibility. In parton terms, 
measuring the small x region at different Q 2 will, 
through analysis of the RG scaling behaviour, en- 
able the gluon distribution Ag(x,Q 2 ) to be iso- 
lated |l2|,^8| accurately. These will be important 
tasks for future experiments at HERA with a po- 
larised p beam. 

7. Semi-inclusive DIS 

As well as pushing further into the small x fron- 
tier, an increasingly important role in the study 
of polarisation and the proton structure will be 
played in future by semi-inclusive processes. We 
have already mentioned the importance of open 
charm production in isolating the gluon distribu- 
tion Ag(x, Q 2 )- In this section, we describe what 
may be learnt from semi-inclusive DIS, both in 
the current and target fragmentation regions. 




Fig. 9 Semi-inclusive DIS: current fragmentation region. 

e 




Fig. 10 Semi-inclusive DIS: target fragmentation region. 

The two distinct contributions to the semi- 
inclusive DIS reaction eN — > ehA from the cur- 
rent and target fragmentation regions are shown 
in Figs. 9 and 10. The current fragmentation 
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events are described by parton fragmentation 
functions Di(jz^, Q 2 ), where i denotes the par- 
ton, while the target fragmentation events are de- 
scribed by fracture functions [p9| M^ N (x, z,Q 2 ) 
representing the joint probability distribution for 
producing a parton with momentum fraction x 
and a detected hadron h (with momentum p' 2 ) 
carrying energy fraction z from a nucleon N. 

The lowest order cross section for polarised 
semi-inclusive DIS is: 

dAa Y P 47ra 2 ^ 
} 1~ 2-^ 



dxdydz 
1 



l-x 



A qi (x,Q 2 )D?(-^,Q 
V 1 — x 



AM™(x,z,Q 2 ) 
(32) 



where z = Eh/ En p0|]. The notation is slightly 
different from sect. 3. Here, Aqi(x) refers to 
quarks and antiquarks separately and a sum over 
both is implied. AM^ N (x, z, Q 2 ) is the polari- 
sation asymmetry of the fracture function. The 
NLO corrections to (32) are given in ref . ifscf . 

Quite different physics emerges from the two 
regions. A programme of semi-inclusive DIS in 
the current fragmentation region has already been 
carried out by SMC 1 31 and will be pursued by 
HERMES and COMPASS. The particular inter- 
est is that it allows the distributions for quarks 
and antiquarks to be separated, giving informa- 
tion on the polarised 'valence' and 'sea' distribu- 
tions defined as Aq^(x) = Aqi(x) — Aq~i(x) and 
Ag| ea (a;) = 2Aq~i(x). These are found by com- 
paring cross section asymmetries for positive and 
negative charged hadrons h, with the assumption 
that the fragmentation functions D\ in (32) sat- 
isfy isospin and charge conjugation symmetry and 
are helicity independent fl32[ . The first results for 
the moments Au 11 and Ad 11 are given in ref. ]3l[| . 

Recently, a new proposal to exploit semi- 
inclusive DIS in the target fragmentation region 
to elucidate the 'proton spin' effect has been pre- 
sented^. The idea is to test the 'target indepen- 
dence' conjecture suggested by the CPV method 
by using semi-inclusive DIS in effect to make mea- 
surements of the polarised structure functions of 
other hadronic targets besides the proton and 
neutron. 

The basic conjecture of ref. j| is that for any 
hadron, the singlet axial charge in (3) can be sub- 



stituted by its OZI value multiplied by a univer- 
sal (target-independent) suppression factor s(Q 2 ) 
determined, up to radiative corrections, by the 
anomalous suppression of the first moment of the 
topological susceptibility yY(0). For example, 
for a hadron containing only u and d quarks, the 
OZI relation is simply a° — a 8 , so we predict: 



12 1 



1 -NS ( 3 1 



(Wi(l + 4.s) a 8 ) 



where 



s(Q 2 ) 



C?(a s ) a°(Q 2 ) 



(33) 



(34) 



Since s is target independent, we can use the 
value measured for the proton to deduce T\ for 
any other hadron simply from the flavour non- 
singlet axial charges. From our spectral sum 
rule estimate of \/x'(0), we find s ~ 0.66 at 
Q 2 = lOGeV 2 , while the central value of the SMC 
result (26) gives s ~ 0.55. 

The non-singlet axial charges for a hadron B 
are given by the matrix elements of the flavour 
octet axial currents, so can be factorised into 
products of SU{2>) Clebsch- Gordon coefficients 
times reduced matrix elements. Together with 
the target independence conjecture, this allows 
predictions to be made for ratios of the first mo- 
ments of the polarised structure functions gf for 
different B. Some of the most intriguing are: 



r?/r? 



2s - 1 - 3(2s + l)i7D 



7r 



2s + 2 - 6sF/D 

if" /if" = 



2s 



2s- 1 



(35) 



(36) 



where (£°) is the state with valence quarks 

uuc (ddc). The results for A and E c are partic- 
ularly striking because of the 2s — 1 denomina- 
tor factor, which is very small for the range of s 
favoured by experiment. These examples there- 
fore show spectacular deviations from the valence 
quark counting (OZI) expectations, which would 
give the ratio 4. 

The proposal of ref.|f| is that these ratios can 
be realised in semi-inclusive DIS in a kinematical 
region where the detected hadron h (a pion or D 
meson in these examples) carries a large target 
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energy fraction, i.e. z approaching 1. While this 
process is most rigorously described in terms of 
fracture functions, it can be pictured as the single 
Regge exchange shown in Fig. 11. The fracture 
function description, which utilises the recently 
introduced extended fracture functions of ref. [ 03[ 
p4j , may be found in ref.||. 

_ e (Pi) 

e (Pi) 




P (P2) 



Fig. 11 Single Reggeon exchange model of ep — > eh A. 

The moment of the polarised structure function 
gf of the Reggeon B can be extracted from the 
polarisation asymmetry of the differential cross 
section moment in the limit z — *■ 1: 

1-2 , dAa tar s et 
dx x — — - — - — — = 

o dxdydzdt 

^^-A/(z,t) / dx B g®(x B ,t;Q 2 )(37) 



z = 
The 



* * Jo 

where z = p' 2 .q/p2-q, x B = Q 2 /2k.q, 1 
x/x B , and t = -k 2 <C 0{Q 2 ) so that z ~ 
emission factor A/(z, t) cancels in the ratios. 

Since the predictions (36) depend only on the 
SU(3) properties of B, together with target in- 
dependence, they will hold equally well when B 
is interpreted as a Reggeon rather than a pure 
hadron state. The ratios (36) can therefore be 
found (see e.g. Fig. 12) by considering the pro- 
cesses ep — > eir~ (D~)X and en — > e7r + (D°)X. 




Fig. 12 Quark diagram for the NhB vertex in the reaction 
ep — > e7T — X where B has the quantum nos. of A^^. 

Of course, the ratios (36) are only obtained in 
the limit as z approaches 1, where the reaction 



cN — > ehX is dominated by the process in which 
most of the target energy is carried through into 
the final state h by a single quark. At the opposite 
extreme, for z approaching 0, the detected hadron 
carries only a small fraction of the target nucleon 
energy and has no special status compared to the 
other inclusive hadrons X. In this limit, the ratio 
of cross section moments for ep — * C7r~(D _ )X 
and en — * e7r + (D°)X is simply the ratio of the 
structure function moments for the proton and 
neutron (35). 

Interpolating between these limits, we expect 
the ratios of f*~' dx x d ^yZdl in the ran § e < 
z < 1 for en — > e7r + (D°)X over ep — > e7r _ (D _ )X 
to look like the sketch in Fig. 13. 






-0.10 
-0.13 
-0.20 
-0.28 
-0.30 



Fig. 13 Cross section ratios for en — > &ir~^~ (D®)X over 
ep — > e7T - (D - )X between Z — > and Z — ► 1, con- 
trasting the OZI and CPV predictions 

The difference between the OZI (or valence 
quark model) expectations and these predictions 
based on our target-independent interpretation of 
the 'proton spin' data is therefore quite dramatic, 
and should give a clear experimental signal. 

Since the proposed experiment requires particle 
identification in the target fragmentation region, 
it is difficult to do at a polarised fixed-target ex- 
periment such as COMPASS |3^], which is better 
suited to studying semi-inclusive processes in the 
current fragmentation region. A better option is 
a polarised ep collider, such as HERA|36[. Test- 
ing our predictions requires comparision of proton 
and neutron data, which can be extracted from 
experiments with polarised deuterons replacing 
the protons in the collider. 
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